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ABSTRACT 

Current theoretical models for the outflows /jets from AGN, X-ray binaries and young stellar objects 
involve large-scale magnetic fields threading an underlying accretion disk. We suggest that such a disk 
is subjected to warping instability and retrograde precession driven by magnetic torques associated with 
the outflow. The growth timescale for the disk warp and the precession period are of order the radial 
infall time of the disk. These effects may be relevant to jet precession and other variabilities observed in 
many systems. 

Subject headings: accretion, accretion disks - MHD - instabilities - magnetic fields - jets and outflow - 
binaries: close 
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1. INTRODUCTION 

Collimated outflows /jets are produced from active 
galactic nuclei, compact stars in X-ray binaries, and young 
stellar objects (YSOs) (e.g., Krolik 1999; Mirabel & Ro- 
driguez 1999; Fender 2003; Reipurth & Bally 2001). Over 
the last two decades, evidence for jet precession in these 
systems has steadily increased. The clearest example is 
the galactic source SS 433, whose jet direction varies with 
an approximate 164 day period (Margon 1984; Eikenberry 
et al. 2001). The black hole binary GRO J1655-40 shows 
jet precession with a period of 3 days (Tingay et al. 1995; 
Hjellming & Rupen 1995). The jets from the super soft 
source CAL83 may be precessing with a period of ~ 69 
days (Cowley et al. 1998). In AGNs, indirect evidence 
for jet precession is seen in the morphology of the radio 
hot spots, which show "fossil" components offset from the 
"present" component positions (e.g., Cygnus A). Finally, 
the changes in the flow directions of several YSO jets have 
been interpreted in terms of jet precession (e.g. Terqucm 
et al. 1999; Bates et al. 2000). 

A natural cause for jet precession is the changes of ori- 
entation of the underlying accretion disk. In addition, the 
super-orbital variabilities observed in a number of X-ray 
binaries (e.g., the 35-day periodicity in Her X-l; see Pried- 
horsky & Holt 1987; Ogilvie & Dubus 2001) have long been 
interpreted as due to precession of a tilted accretion disk. 
In both X-ray binaries and YSOs, the tidal force from the 
companion star could play a role in driving the preces- 
sion (e.g., Katz 1973; Wijers & Pringle 1999; Terquem et 
al. 1999; Bate et al. 2000; Ogilvie & Dubus 2001), although 
it may not be the only or the dominant effect (e.g., the 3 
day precession period of the GRO J1655-40 jet is too short 
to be explained by the tidal effect for a 2.6 day binary). 
Most importantly, for the precession to operate, disk tilt 
needs to be excited and maintained. In accreting binary 
systems, the disk plane is expected to be aligned with the 
orbital plane since the disk angular momentum originates 
from the binary motion. For YSOs in binaries, the disk 
plane may be initially misaligned with the orbital plane. 
However, if we consider the disk as a collection of circular 
rings, different rings will have different precession rates; 
it has been recognized that the combined effects of differ- 



ential precession and internal disk stress/dissipation tend 
to damp the the disk tilt, so that the disk settles into the 
binary plane (Lubow & Ogilvie 2000). 

Several driving mechanisms for disk tilt /warp have been 
proposed. Schandl & Meyer (1994) showed (in the con- 
text of Her X-l) that irradiation-driven wind from the 
outer parts of the disk can carry away momentum flux 
and make the disk unstable to warping. Pringle (1996) 
showed that even without wind loss, radiation pressure it- 
self can induce warping instability in the outer region of 
the disk. Futher studies indicated this radiation-driven 
warping instability may indeed operate in X-ray binaries 
(e.g., Maloney, Begelman & Nowak 1998; Wijers & Pringle 

1999) , although it does not provide a generic explana- 
tion for the long-term variabilities in all X-ray binaries 
(Ogilvie & Dubus 2001). Quillen (2001) showed that a 
wind passing over the disk surface may induce warping via 
Kelvin-Helmholtz instability. Finally, in the case of disk 
accretion onto magnetic stars (e.g., neutron stars, white 
dwarfs and T Tauri stars), the stellar magnetic field can 
induce disk warping and precession (Lai 1999; see also Aly 
1980; Lipunov & Shakura 1980; Terquem & Papaloizou 

2000) ; this may explain several observed features of quasi- 
periodic oscillations in low-mass X-ray binaries (Shirakawa 
& Lai 2002a), milli- Hertz variabilities in accreting X-ray 
pulsars (Shirakawa & Lai 2002b), and variabilities of T 
Tauri stars (Terquem & Papaloizou 2000; see also Agapi- 
tou et al. 1997). 

In this paper we suggest a new disk warping mecha- 
nism that is directly tied to the production of magnetically 
driven outflows/jets. Using an idealized setup (§2), we 
show that a disk threaded by a large-scale magnetic field 
may experience a warping instability and precess around 
the central object (§3). These magnetically driven disk 
warping and precession arise from the interaction between 
the large-scale magnetic field and the induced electric cur- 
rent in the disk. While more studies are needed, we sug- 
gest that these effects may provide a natural explanation 
for the procession of jets/outflows and other variabilities 
observed in various systems (§4). 

2. THE SETUP 
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The current paradigm for the origin of astrophysical 
jets/outflows involves a large-scale magnetic field thread- 
ing the accretion disk around a central object (star or 
black hole); this ordered magnetic field plays a crucial role 
in extracting/channeling mass, energy and angular mo- 
mentum from the disk. The energy outflow can be either 
hydromagnetic (with significant mass flux) (Blandford & 
Payne 1982) or electromagnetic (dominated by Poynting 
flux) (Blandford 1976; Lovelace 1976). The origin of the 
disk-threading magnetic field is not completely clear: the 
field could be advected inwards by accretion, or generated 
locally by dynamo processes. In the case of protostellar 
outflows, the stellar magnetic field may play an important 
role (Shu et al. 1994, 2000). Many theoretical/numerical 
studies have been devoted to understanding magnetically 
driven outflows/jets from accretion disks (e.g., see recent 
reviews by Lovelace et al. 1999; Konigl & Pudritz 2000; 
Meier et al. 2001). 




Fig. 1. — A sketch of the idealized disk - magnetic field configu- 
ration in which hydromagnetic outflows are produced. The disk is 
threaded by a large-scale poloidal field, which gets twisted by the 
disk, generating the toroidal field. 

Figure 1 shows an idealized sketch of the magnetic field 
- disk configuration in which hydromagnetic outflows are 
produced. A geometrically thin disk is threaded by the 
poloidal magnetic field B p = BrR+BzZ (where R, Z are 
unit vectors in cylindrical coordinates). This poloidal field 
is of course curved (with Br = at the disk midplane), 
and we use the superscript "+" to denote the field at the 
upper disk surface. For a sufficiently inclined field (ra- 
dial pitch angle = tan -1 \B\jBz\ > 30°), centrifugally 
driven outflows are possible (Blandford & Payne 1982). 
The twisting of B p by the disk and the outflow gives rise 
to a toroidal field which has different signs above and 
below the disk plane. We introduce the azimuthal pitch 
C, so that B^, = —(Bz at the upper disk surface. That 
( is positive reflects the fact that the disk always rotates 
faster than the outflow 1 . The discontinuities in £?$ and Br 
across the disk imply a net disk surface current (integrated 
over the disk thickness) K = (c/2tt)(- j B+R + B+&). 
Thus, for an unpertutbed (unwarped) disk, the only torque 
on the disk (per unit area) is the usual magnetic braking 
torque 

T z = ^-RB z B+ = -^-(RB z . (1) 



2tt"~" * 2tt 

1 In the stationary, axisymmetric, MHD wind theory, B$/B p = R(Cl 
and Q(R, Z) its rotation rate, fla is the disk rotation at the midplane. 
schematically by the equation dBg>/dt = RBzdQ/dZ — -B^/rdiss, where 
with £ = (R/HB^Tdisg (where Hb is the vertical scale height in which 



Angular momentum is extracted out of the disk by the 
open field lines, and can be carried away by torsional 
Alfven waves or a centrifugally driven wind (for sufficiently 
inclined poloidal fields). In addition, the radial magnetic 
force, 

F R = B z B+/(2n), (2) 

provides partial support of the disk against gravity. 

The configuration shown in Fig. 1 lies at the heart of 
many theoretical models of magnetic disk outflows. The 
real situation is certainly much more complicated. Nev- 
ertheless, in §3 we shall adopt this simple picture in our 
consideration of magnetically driven disk warping and pre- 
cession — our goal is to discuss these physical effects in 
a way that is as transparent as possible, so as to expose 
their robustness and potential pitfall. 

3. MAGNETICALLY DRIVEN DISK WARPING AND 
PRECESSION 

We now consider warping perturbations to the disk. The 
disk can be considered as a collection of circular rings 
which interact with each other via internal stresses. We 
specify the disk warping by the unit normal vector l(r) for 
the ring at radius r. Consider the situation (see Fig. 2) in 
which the ring around radius r is tilted by an angle while 
most of the disk remains flat (li = I2 = Z, where Z is the 
normal vector of the unperturbed disk; see Fig. 1). Since 
the magnetic field Bz is produced by all the currents in the 
disk-outflow system, we posit that Bz at radius r is not 
modified by such a local perturbation. For the ring at r, 
the field Bz acts as an "external" magnetic field. Now set 
up a coordinate system with the z-axis along the l(r) (see 
the lower panel of Fig. 2) and the associated cylindrical co- 
ordinate system (r<j>z). As in the case of the unperturbed 
disk, twisting of B z = B z cos by the disk rotation pro- 
duces a toroidal field (in the new coordinate system) B^^ 
(where the subscript "t" serves as a reminder that this 
field is due to twisting of B z ), whose value (at radius r) 
at the upper disk plane is Bf t = —C,B Z) and at the lower 



disk plane B , 



(,B Z . Note that B z is necessarily contin- 



uous across the disk plane in this perturbation. The radial 
surface current associated with (B~f t — B~7 t ) is 



K r 



"277^* 



2tt 



C,B Z cos 0. 



(3) 



Interaction of this K r with B z gives an azimuthal force 
on the disk, and results in the magnetic braking torque 
(per unit area) T z = -(rB 2 /(2ir) = -(rB 2 z cos 2 0/{2ir), 
analogous to cq. (1). Unlike the unperturbed disk, how- 
ever, the "cxtrcnal" field H z now has a ^-component, 
B(f, = B z sin0 cos<j), and the interaction of K r with B$ 
gives rise to a vertical force (per unit area) 



F z 



1 

2^ 



C B\ sin cos cos ( 



(4) 



This vertical force can also be seen as arising from the 
imbalance of magnetic stresses above and below the disk 
plane: B^ = B^ + B~£ t in the upper disk plane while 

— Q.d)/V p , where V P (R, Z) is the the poloidal velocity of the outflow 
Alternatively, in the absence of outflow, B$ near the disk is governed 
Tdiss is the disspation timescale. In steady state, we find B^ = —C,B Z , 
fi in the magnctosphcre varies. 
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£>0 = Bj, — B^ t in the lower disk plane, and thus 

F z = [(B^) 2 -(B^) 2 ]/(8tt). Obviously, averaging the force 

over the azimuthal angle, (F z ) = [1/2-k) L" F z d(j>, gives 
zero. But since F z is unevenly distributed along the ring, 
it gives rise to a net torque on the ring. Averaging over <p, 
the torque per unit area is 

(T warp ) = --!-£ r B 2 Z cos sin y 



4ir 

47T 



C rB% cos 



(z-i) 



(5) 



Clearly, this torque tends to pull the vector l(r) away from 

"external" magnetic field direction Bz. This is the mag- 
netically driven warping instability: for small tilt, the angle 
tends to grow (in the absence of other forces) at the rate 
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where £ is the surface density, and fid the angular veloc- 
ity of the disk. In a hypothetical situation, if had a 
fixed direction for all 0, and if there were no coupling be- 
tween different rings, the vector l(r) would evolve toward 
being perpendicular to B^, i.e., the plane of the disk (ring) 
would prefer to lie along the "external" magnetic field (cf. 
Lai 1999). 





Fig. 2. — Warping perturbation to the disk around radius r. For 
clarity, only the Z-component of the magnetic field, B^, is shown. 

In the perturbed disk, the azimuthal surface current 
K$ ~ = (c/2ir)B+ = (c/2ir)B z tan0 (for small 0) 
also interacts with the "external" field B^, which can 
be decomposed into B^ = Bz cos z + Bz sin (sin f + 
cos 4> 4>) ■ The interaction between K$ and B z = Bz cos 
gives a radial force, as in cq. (2). For the perturbed disk 
(/3 / 0), a new vertical force arises from the interaction 
between Ka, and B r = Bz sin sin <j>: 



This force is also equal to [(B~) 2 - (B+) 2 ]/(87r), where 
Bf = B z sin sin (f> ± B z taxi6. Again, (F z ) = 0, but 
the net torque is non-zero. Averaging over <j>, we find the 
torque per unit area 

(T prcc ) = ~K+B x r sin/3x = ~ K^B z rt x 1. (8) 

Clearly, the effect of this torque is to make the ring (at 
radius r) precess around the Z-axis (the unperturbed disk 
normal vector). The angular frequnecy of the magnetically 
driven precession is 

_ K,B Z ± B 2 tzne ± 

The negative sign implies that the pecession is retrograde, 
i.e., opposite to the disk rotation. 

4. DISCUSSION 

To get an idea on the warping and precession timescales 
associated with eqs. (6) and (9), we note the disk angu- 
lar momentum equation (for an unwarped disks) can be 
written as 



(10) 



where Vr is the radial velocity and T v ; sc is the viscous 
torque. If we set T v j sc = (/ — l)(—(RB z /2ir) and use 
n d ~ (GM/R 3 ) 1 / 2 , we find 



warp 



C in \ = - Vr 



tan 6 



(11) 



■ K^Bz sin/3 sin< 



(7) 



AfR 

For typical parameters (C ~ tan^ ~ / ~ 1), the warping 
timcscale and precession period are of order the disk infall 
time R/\V R \. 

It is important to note that the magnetically driven disk 
warping and precession discussed in §3 are secular effects 
which operate on timescales much longer than the disk 
dynamical time. In deriving eqs. (6) and (9), we have 
implicitly assumed that the perturbed disk (see Fig. 2) 
reaches a new steady state so that the magnetic field be- 
haves in a similar way as in the unperturbed disk (e.g., 
B^ t = =f££? 2 ). These effects are not captured in a dy- 
namical perturbation analysis of the disk. Clearly, our 
derivation of the magnetically driven warping and preces- 
sion is far from being rigorous. Nevertheless, it suggests 
that the warped, processing disk-outflow may be an alter- 
native (and preferred) state for the accretion disk threaded 
by large-scale magnetic fields. 

The effects of magnetically driven warping and preces- 
sion discussed in this paper are analogous to the similar 
effects that exist in the accretion disks around magnetic 
stars (Lai 1999; Shirakawa & Lai 2002a, b; Pfeiffcr & Lai 
2003; see also Aly 1980; Lipunov & Shakura 1980). In 
those systems, the central star provides the "external" 
magnetic field for the disk. The warping and precessional 
torques arise from the interaction between the surface cur- 
rent on the disk and the horizontal magnetic field (parallel 
to the disk) produced by the stellar magnetic dipole. Here, 
in the absence of the central magnetic star, the "external" 
magnetic field (for the warped region of the disk) is gen- 
erated by the other currents in the disk-outflow system, 
and we have argued that similar torques exist: the warp- 
ing torque relies on the surface current generated by the 
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twisting of the vertical field threading the disk, while the 
precessional torque relies on the azimuthal screening cur- 
rent due to the diamagnctic response of the disk. The 
simplistic nature of our arguments (perhaps even to the 
point of being wrong) in §3 and the complexity of the real 
astrophysical systems (see references in §2) preclude us 
from addressing many important issues such as the extent 
and location of the warped disk region. Nevertheless, we 
are not aware of any previous discussion of these effects, 
and we think it is useful to present them so that they may 
be studied further by the astrophysics community. 

The magnetically driven warping and precession effects 
discussed in this paper potentially have advantages over 
the other disk warping mechanisms that are needed to ex- 
plain jet precession (see §1). For example, the radiation 
driven warping instability (Pringle 1996), while probably 



important for X-ray binaries (Ogilvie & Dubus 2001), is 
irrelevant for YSOs (since the instability operates at large 
disk radii). The magnetic effects discussed in this paper 
are directly connected to the large-scale magnetic field as- 
sociated with the jet/outflow, and the model naturally pre- 
dicts retrograde precession of the warped disk (and thus 
the jet), in agreement with observations in X-ray binaries 
(e.g., Her X-l). 

The idea for this paper was motivated by the questions 
from Chris McKee and Frank Shu following a seminar that 
the author gave at UC Berkeley in 2001. I thank Richard 
Lovelace for useful discussion. This work is supported in 
part by NSF Grant AST 9986740 and NASA grant NAG 
5-12034, as well as by a research fellowship from the Alfred 
P. Sloan foundation. 
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